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Theory of massless scalar field with interaction g<p 3 in six-dimensional space 
is considered. A possibility of initial scale invariance breaking, which results in 
a spontaneous arising of effective interaction G<p A 1 is studied by application of 
Bogolubov quasi-averages approach. It is shown, that compensation equation for 
form-factor of this interaction in approximation up to the third order in G has 
a non-trivial solution. In the same approximation Bethe-Salpeter equation for a 
zero-mass bound state of two scalar fields is shown to have a solution. The 
conditions imposed on form-factor value at zero and scalar field mass m fix the 
unique solution, which gives relations between parameters of interaction g<p 3 and 
parameters G and m. Arguments are laid down in favour of a stability of the 
non-trivial solution. 
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1 Compensation equation in quasi-averages approach 

N.N. Bogolubov quasi-averages method |2j is the most consistent and effective method of 
studying of a spontaneous symmetry breaking problems. 

An important point of the quasi-averages method is connected with a compensation 
equation |2] . Bearing in mind applications in the present work let us briefly formulate 
method of construction of the compensation equations. In the line of a study of a possible 
spontaneous symmetry breaking in quantum field theory problems in method [2] the following 
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procedure is applied 1 . Let the initial Lagrangian 

L = L + L int ; (1) 

to possess some symmetry. Let us add to expression (0) some term eL^, which breaks the 
initial symmetry. With this modification of the problem we perform evaluations of necessary 
quantities and we set e — > only after these evaluations. Not always the results of such 
a procedure (quasi-averages) coincide with results, obtained in the framework of the initial 
symmetric problem (simply averages). In the line of these evaluations of quasi-averages 
one has to solve compensation equations. For instance, in a theory with the initial chiral 
symmetry fermions are to have zero masses. Let us use the following small increment which 
breaks the symmetry 

eL br = - • ( 2 ) 

Now let us add to the modified Lagrangian (J2J) a possible mass term and subtract the same. 
We have 

L = L - rmjjijj + L int + m-0-0 - e^-0 (3) 

Let the fist two terms to be the new free Lagrangian while the three last terms now comprise 
the new interaction Lagrangian. Then we have to demand the new interaction does not 
contribute to the mass term, that is two-field Green function obtained from the modified 
interaction Lagrangian be zero on the mass shell. This condition is just the compensation 
equation of the problem. In the case under consideration this condition leads to equation 

- m + e + S(m) = 0; (4) 

where S(m) is mass operator on the mass shell of the modified free Lagrangian. In this 
equation one already can set e — > 0. As a rule (see e.g. [3]) mass operator E(m) is proportional 
to m and trivial solution of the compensation equation m = always exists. However 
nontrivial solutions m ^ also may exist. 

Thus the main principle of construction of a compensation equation consist in the pro- 
cedure "add - subtract" of symmetry breaking terms, one of these terms being related to 
the free Lagrangian and the other one being related to the interaction Lagrangian. Then 
one has to compensate that term, which is to be zero in the corresponding problem. This 
principle will be applied in the present work. 

2 Justification of the model choice 

The phenomenon of spontaneous symmetry breaking is decisive for formulation of the elec- 
troweak theory. The introduction of elementary scalar fields j3] is well-known to be the 
essence of the standard way of the breaking. Their self-interaction leads to redefinition of 

1 At first methods where applied to quantum theory problems in work [3] 
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the vacuum in the theory and to existence of the scalar Higgs particle. However proposals 
based on a dynamical breaking of the initial symmetry without elementary Higgs scalars are 
also considered. This leads to effective four-fermion interaction of heavy quarks either to 
be postulated (see e.g. review [5j) or to be dynamically arisen [HI El- As a result the initial 
theory, which contains scale-invariant gauge interactions, transforms into a theory, which 
contains interactions with dimensional coupling constant as well, that explicitly breaks the 
scale invariance. In works being mentioned |H1 U\ assumption was made on a possibility of 
existence of solutions of corresponding compensation equations. However this assumption 
was not duly justified. The purpose of the present work is to consider a simpler model, 
which would allow to have exact solutions of (approximate) compensation equations. Using 
these solutions one could study conditions under which the assumptions would be fulfilled. 
To some extent the model has to correspond to features of the electroweak theory. Namely 
we achieve a simplicity by considering a scalar field. In view of coupling constants to have 
proper dimensions we choose dimensionality of the space-time to be six. Really, in this case 
the coupling constant of interaction g (ft 3 is dimensionless and interaction G (ft 4 has constant 
of inverse mass squared dimension, that corresponds to the dimension of a constant of a 
four-fermion interaction in four- dimensional space-time. 

So we introduce in the six-dimensional space-time a scalar field (ft with initiative scale- 
invariant Lagrangian 

1 d(j> d(ft g 3 
L= 2 9 d^ + 3! - (5) 
Let us choose the natural signature with one time and five space axes. The transition 
from this space-time to Euclidean six-dimensional space is accompanied by the following 
substitutions 

p2_>_p2. d 6 P ^ld% P . (6) 

It was important for us to find a model, which corresponds to the approach under consid- 
eration. So in the present work we will not discuss physical meaning of a multi-dimensional 
theory and we consider the chosen variant as purely model one, as well as two-dimensional 
models are often considered. 

Now we start with Lagrangian (|5|). Evident evaluations give one- loop renormalization 
group equation [8 J for g 2 (fi 2 ) 

dL 4(4tt) 3 ' g A 3 



£ = log^; (7) 



Solution of equation J7J) has a form 



2/ 2\ 2 ( -i ^ g 2 n 
9^)=9o{^ + JU-TE log 



4(4vr) 3 °A 



Sometimes it is convenient to use parameter h(fi 2 ) defined by the following relation 

4(4tt) 3 V 



W> = S = (><*£) ; (9) 
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where for transition from A3 to A 2 , we have used the standard tool analogous to that in QCD: 



Thus we get convinced, that the theory (J3J) is an asymptotically free one and expression (J£JJ) 
makes sense for /z 2 ^> A 2 ,. 

Note that in this theory there are quadratic divergences in the scalar field mass. It is the 
common feature of theories with elementary scalars. The problem of mass of the scalar field 
will be considered in details later on. 

3 Compensation equation in a six-dimensional scalar 
model 

Let us have a massless scalar field of the six- dimensional space. The initial free Lagrangian 
possesses scale symmetry. We shall look for a solution, which breaks this symmetry, with 
the aid of Bogolubov quasi-averages approach 0. Namely according to [2] we add to the 
Lagrangian the following small increment 

Now the scale invariance is already broken and an appearance of nonlocal terms of the form 

G J F(x, xi, Xz, x 3 , £4) </>(xi)<f>(x2)</>(xs)^)(xi) dx\ cfo 2 dx 3 dx^; (10) 

is possible. Here G is a dimensional coupling constant and F(x, x\, X2, £3, £4) is a function 
of four differences of coordinates x — X{, which Fourier transform F(j>i,p2,pz,Pi), where pi 
are momenta of legs, represents a form-factor, defining range of interaction (fTU|) . We shall 
look for a solution, decreasing at momentum infinity and thus defining a region of action of 
the effective interaction. 

Let us add to the initial free Lagrangian such a term with an interaction of the forth 
power and subtract the same 

2 y dx^ dx v 2 v 4! v 4L 3! v 4! v 2 v ' v ; 

where we use abbreviated notation — G F ■ (ft 4 instead of expression (|10p. Of course the 
presence of term (jlUJl explicitly breaks the scale invariance, so we perform a procedure " add 
- subtract" for a mass term as well. Let us refer the forth power term with the plus sign 
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to the interaction Lagrangian and the same term with the minus sign we refer to the free 
Lagrangian. 

1 d(f) d(j) m 2 2 G 4 e 4 
L mt = |0 3 + ^F.0 4 + ^0 2 ; (12) 



According to the quasi-averages approach [2 j the term with the plus sign has to be com- 
pensated. This means, that the new free Lagrangian leads to zero four-particle connected 
Green functions and as a final result contains only terms of the second power in fields. Thus 
performing evaluations with sign which is inherent to the term in the new free Lagrangian, 
we come to the compensation equation, which schematically looks in the following way: 
the first order term plus one-loop terms plus two-loop terms etc. Emphasize once more, 
that here one has to use term + G 4 as an interaction Lagrangian. One has to equalize 
to zero the expansion obtained in such a way . This condition is an equation for function 
F(pi,P2,P3,P '4)- We set e — > after evaluations, in our case this means after compensation 
equations being obtained. 

The equation explicitly differs from expansion in powers of interaction Lagrangian 

Unt = ^F-^- (13) 

in the sign of the interaction constant. In view of this note let us emphasize, that the 
procedure being described can be applied only to symmetry breaking terms of even powers 
in fields. For terms of odd powers, e.g. for three-linear ones, a fulfillment of a compensation 
equation leads to vanishing of connected Green function, which is defined by an interaction 
Lagrangian, because the two expansions in this case differ only in overall sign. 

Note, that the presence of term — G 4 in the new free Lagrangian may lead to appearance 
of connected Green functions of higher powers in 0, that is of the sixth power, of the eighth 
power etc. Generally speaking, one has to construct a chain set of compensation equations 
for all these Green functions. We start with an equation for the fourth power Green function 
and the problem of higher Green functions will be discussed in what follows. 

Let us construct an approximated equation for the fourth power connected Green func- 
tion. First of all we choose the following kinematics: both left legs have zero momenta and 
the right ones have momenta p and —p. We restrict ourselves by terms up to two-loop ones 
inclusively. Namely, we have the first order term - the point; three terms of the second 
order - simple loops, i.e. a horizontal one and two vertical ones with permuted left legs; in 
the third order we have a horizontal and two vertical two-loop chains and six terms "wine 
glass": horizontal wine glasses having bases to the left and to the right and vertical ones 
with bases up and down. The number of the last terms is to be counted twice due to permu- 
tations of the left-sided momenta p and —p. Generally speaking, in each vertex form-factor 
F is present. However we can solve only a linear version of the equation, which is obtained 
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by keeping in the equation the first and the second order terms, the two-loop horizontal 
chain and the wineglass with the basing to the right. Contributions of the rest third or- 
der terms we shall consider later on. We proceed to the linear equation keeping form-factor 
F(p, —p, 0, 0) = F(p 2 ) in the first order term and in right-hand vertices of the horizontal loop 
of the second order, of the horizontal two-loop chain and of the wineglass in the third order. 
Other vertices in diagrams we consider to correspond to point-like interaction in which the 
form-factor is changed for its value at zero (F(0) = 1) 

§W = §^. (14) 

In vertical simple loops, which as well serve as a kernel of the integral equation, we substitute 
point-like vertices (fPfj) . Corresponding integrals diverge of course. In view of our search for 
decreasing solutions at momentum infinity for F(p 2 ), we introduce some cut-off A, which 
existence is to be confirmed by results of a solution of the equation. In doing this we make 
the following substitution 

pOC 

dq 2 - A 2 . 



o 



For estimation of A order of magnitude we use the following definition 

poo 

A 2 = / F(y)dy; (15) 



where one of vertices is changed for the form-factor. For justification of the approach the 
problem of convergence of the integral in ()15|) . We shall use the same cut-off A in logarith- 
mically diverging integrals. A possible difference of an actual cut-off in these integrals from 
A leads to some change in constant term cj which enters into corresponding expressions. It 
will come clear, that the solution will not depend on a value of this constant. Thus the for- 
mulation of the equation in the framework of the accepted approximations does not contain 
arbitrary assumptions. 

We consider the equation in six-dimensional Euclidean space with the aid of substitu- 
tions (jHJ). In the course of evaluations one has to perform angle integrations in six- dimensional 
space of functions ((p — q) 2 )^ 1 and log(p — q) 2 with powers of (pq). We have (for the loga- 
rithmic case see [§]) 

/ p 2 + q 2 -2pqcos9 = IT ~~ V) (i^ ~~ 4^) + ^ ~ %) (% ~~ v)) ; 

J dtt 6 \og(p 2 + q 2 -2pqcos6) = ^ ^6(x - y) - ^ + 12 log x) + 

+ e( 1 /-a:)(— -4 + 12 1ogy)V ( 16 ) 



/7r 3 / /3y 2 
dtt 6 (pq) \og(p 2 + q 2 - 2pq cos6>) = — I Q(x - y) - 6 



3y 2 _ 3y*\ 
— * y -^) + 
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,/3x 2 „ 3a; 3 \\ 



X = p ; y = q ■ 

First of all let us calculate one- loop integral keeping terms of zero and the first orders in m 2 . 
We have for one such vertical diagram (x = p 2 , where p is the total momentum along the 
loop) 

- ( A2 + r log (^) +2m2 log (^) - ") ■ (17) 

where A is the square of the cut-off being mentioned and c is a constant, depending on a 
behaviour of the form-factor. 
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Fig. 1. The graphic representation of the linear compensation equation (18). 



Let us consider the linear compensation equation, obtained in agreement with the for- 
mulated rules (see Fig. 1). The equation in this approximation has the following form 



GF(p 2 ) 
G 3 



G 1 



2(4tt) ; 



(SA 2 + \p 2 log(^) + 4m 2 log(^) - 2cp 



x 



8 (2tt) 9 

F{q 2 ) 
(q 2 + m 2 )'- 



\{P - Q) 2 log + 2 m 2 log ^£ - c (p - q) 2 ) 



x 



d 6 q - 



3C 3 7r 3 A 2 
2(2 



7T A /• 



(g 2 + m 2 ) : 



(18) 



Firstly let us note, that trivial solution G = is evidently possible. In view of looking for 
a non-trivial solution we cancel the equation by G. Performing here angle integrations by 
using formulas we obtain the following one-dimensional integral equation 



F(x) = (3 A 2 + \x log(iL) + 4m 2 log(^) - 2cx ) - 

3G 2 A 2 [°° y 2 F(y) J G 2 f 1 f x y 5 F(y) 



d y - w7-^ - / 7—^ d y + 



4(4tt) 6 J (y + m 2 ) 2 * 18(4tt) 6 V. 20x 2 ,/ (y + m 2 ) 2 

3 f y 4 F(y) [ x Jl£{y) [* y 2 F(y) A 

+ T~ / 7 — I 2^2 "2/ + 3 log a; / - — ■ — + 3 x log x / - — ■ — dy + 

4^ Jo (y + m 2 ) 2 J (y + m 2 ) 2 J (y + m 2 ) 2 

f°°f]ogyF(y) f°° (4 + 3 log y)y 2 F(y) ^ 
+ 3 / -7 — ■ dy + x ? — ■ ^- dy + 

Jx (y + m 2 ) 2 J x (y + m 2 ) 2 

A C y 3 F(y) , 3x 2 f°° yF(y) , x 3 f°° F(y) \ 
+ 4 / f— Ky ' dy + — / -1-m-dy-— l % dy + 19 

Jo (l/ + m ) 4 A (2/ + ™ 2 ) 2 20 ^ (y + m 2 ) 2 J 

G 2 m 2 ( 1 f y 4 F(y) 8 f y 3 F(y) A* y 2 F(y) 

12(4vr)6^ x 2 7 (y + m 2 ) 2 y+ xJ (y + m 2 ) 2 V + gX J (y + m 2 ) 2 V + 

p y 2 log yF(y) ^ P yF(y) , 2 f°° F(y) , \ 



G 2 

+ 



6(4vr; 



(y + m 2 ) 2 ■ (y + m 2 ) 2 (y + m 

A a 2 1 o \ ( f°° y 3 Hy) A , r y 2 F(y) , \ 

l logA+3c )U ¥T^y dy + x l ^+-m^ 2dy ) + 



G 2 m 2 [°° y 2 F(y) 

(4tt) 6 J (y + m 2 ) 2 

A method of solution of equations of (fl~9"j) type is developed in work [TTjj. Equation (fl"9"j) 
is reduced to a differential one by sequential differentiations. Evident evaluation gives 

(x 2 F(x ))] - p( F ^ I 2m 2 fx ^ F ^ \3 d F ^ )\ 
dx 4 \ dx 4 V J J \(x + m 2 ) 2 V dx 2 (x + m 2 ) 2 dx (x + m 2 ) 2 J J ' 

9 G 2 

" = Wf- (20) 

One easily see, that Eq. (|2U|) can be rewritten in the form 

dx ^ ) ( dx ^ ) ( dx) ( dx) ( dx ) ( dx ) 

xfx A_ 2 N )fx-^--3 > ) + /5x 2 V(^) = 2(3m 2 x(F(x)+x^-x 2 ^- F -); (21) 
V dx / V dx / } V ax ax 2 / 



where two terms of expansion in m 2 are kept. From this form of the equation we immediately 
conclude, that for x — > there are eight independent asymptotes, which coefficients we 
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denote as follows 



a_ 2 a_i 

— ; a ; a ; log x; a x x; 

x^ x 



x log x; a 2 x 2 ; a 3 x 3 . (22) 
Eight independent asymptotes at infinity are the following 

F k (x) ~ x" 3 / 8 exp^ (^x^^exp ^^" 1 ^ ; fc = 1, 2, 8. (23) 

Four of these asymptotes at infinity decrease exponentially (k = 3, 4, 5, 6), and the rest 
four ones do increase. 

Equation (}2~T]) is equivalent to the initial integral equation under definite boundary con- 
ditions being fulfilled. First of all we can use only solutions, decreasing at infinity. To obtain 
conditions at zero we have to substitute expression 

F(x) = -%r4-r( x 2 -^—(x 2 F(x) 



(3 dx A \ dx 4 

in integrals of equation (fHJj) and perform sequential integrations by parts. The results are 
presented in the Appendix. 

Substituting expressions (|57j) into equation ([!§)). we have 

d 1 

Fix) = Fix) — — 7T — — a i log x — aux log x + 

Gtt 3 / t2 / GI 
+ 3A J 1 



2(2vr) 6 V V 2(4 vr) 3 

/ = r y2F{V) dy 
(y + m 2 ) 2 



T l0g (^) - 2cx ) + ^(logA 2 + 3c)a^;(24) 



From here we obtain the following condition (independently on values of A 2 and c) 

2Gm 2 



a_ 2 = 0, a_i = 0, a i 



(4tt) j 



ail = 3(2^ = "a" 5 (25) 

The first four conditions (J25j) are boundary conditions for Eq. (I21|) . A combination of four 
solutions decreasing at infinity with account of these boundary conditions gives the unique 
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solution. It can be expressed in terms of well-known special functions for case m 2 = 0. 
Indeed, let us make the following substitution in Eq. (|2"T]) 

^ = ^r- ( 27 ) 

which reduces the equation to the canonical form of Meijer equation of the eighth order 
d + l)(z± + I)(z4-)(z±)(z4- - l)(z± - 



dz J \ dz 2J \ dzJ \ dz/ \ dz 2 J V dz 2 

"('s-Ohs-D + O^ - (28) 

Conditions (|25jl fix the solution. Firstly, four solutions, decreasing at infinity, always could 
be combined to set to zero three singular asymptotes at zero, i.e. to fulfill conditions a_2 = 
a_i = clqi = 0. Such property has the following Meijer function (see [T2] ) 

C ■ Gll( z | 3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1). 

The constant is defined by the coefficient afore \/ilog z. For small z this Meijer function is 
such [TJ] 

16 _ 

G 5 ° 8 (z | 3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) = vr + — ^ log z + ... . (29) 
Comparing the coefficient afore yfz log z with (|25|1. we obtain 

Performing integration (see (ISj), we have in accordance with definition of J (|24|) 

/OO I7y l>00 r\ fey 

F(y) dy = ^J Q Gl° 8 ((3y 2 /2 S \ 3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) dy = (30) 

that perfectly agrees with condition (|2 
Thus, solution 

F(x) = ^ G™(f3x 2 /2 8 | 3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) ; (31) 

fulfills all conditions (|25|). and consequently the initial equation (|19jh which is an approximate 
compensation equation. This solution is a nontrivial solution, which contains dimensional 
parameter G, and hence it leads to the initial scale symmetry breaking. Of course as we have 
noted before trivial solution F(x) = is also possible. Note, that the boundary conditions 
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are not dependent on value of the form-factor at zero. Equality .F(O) = 1 will serve as an 
additional condition in what follows. 

Let us take into account terms proportional to m 2 . We shall look for a correction to the 
solution of Eq. (|21|) in the following form 

F(x) = F (x) + AF(x) . (32) 

Substituting (|32j) into equation (J21j) we have the following equation in the first order in m 2 

dx / ( dx ) ( dx) ( dx) ( dx ) ( dx ) ( dx ) 

x { x Tx ~ 3 ) + ^) AF{x) = 2 ^ m2x ( Fo{x) + xd -^- x2 ^)- (33) 



From equation ()33j) we can exactly define several terms of expansion of A F(x) for small 
x. Indeed let us consider the following expression 

T . , 2 m 2 / . , d F 9 d 2 F \ , , 

AF(x) = (F (x) + x — - - x 2 — - ) = (34) 



x V dx dx 2 

n i( 71 2 G / . /r . 23\ irV2G 2 , ^. A 

= 2m 2 — =- H , — — log( v / /5x) +47 H — -i log x + O(x) ; 

\2V2x 3(4tt) 3 V BKVh ' 1 6/ 96(4vr) 6 6 v 7' 

where 7 = 0.577215665... is the Euler constant. Substituting expression (jHljl into equa- 
tion we get convinced, that it fulfills the equation up to terms of x 3 order, because 
the differential operator in the left-hand side nullifies the terms presented in ()34j) and sub- 
sequent terms up to the indicated order. We are interested just in the presented terms (l3~H) 
because they refer to the boundary conditions. Indeed expression (J34|) contains terms 
z -1 / 2 , log z, z 1 / 2 log z, which violate their boundary conditions. Hence we are to add to 
expression (|34|) a combination of solutions of the homogeneous equation to force the bound- 
ary conditions to be fulfilled. Finally we obtain 

AF(x) = A F(x) - ^0/3x 2 /2 8 |3/2, 1, 1/2, 1/2, -1/2, 0, 0, -1) - 
-^-Gl° 8 ((3x 2 /2 8 \3/2, 1, 1/2, 0, 0, 1/2, -1/2, -1) - 

- tt Y (7 + log 2 - G™(f3x 2 /2 8 I 3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1 ) ; (35) 

Gm 2 
~ 274^)3' 

From this expression we extract the exact value for F(0). While doing this one has to bear 
in mind, that the presence of a term being proportional to log x at x — >• is a consequence 
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of an expansion in m 2 at x ^> m 2 . Looking back at the corresponding evaluations we see, 
that for x — > one has to change log x for log 4 m 2 . Now we have 

F(0) = + y(41ogr+(16-7r 2 )7 + (14-7r 2 )log2-^ + 7r 2 ^). (36) 

4 V 9 48 / 

For Y = we obtain F(0) = 1.11072. Condition F(0) = 1 defines the value of Y, 
which is connected with mass (see (J3*5|) ) 

Y = 0.005789. (37) 

Thus the solution, which is found here, satisfies all the necessary conditions provided (|3*?j) is 
valid. Emphasize, that (J37|) defines the mass of the scalar field. Note, that the small value 
of (J37|) thoroughly justifies the account of only the first term of the expansion in m 2 . We 
reject the second solution of condition F(0) = 1, which is of order of unity, due to to its 
inconsistence with the expansion of the solution in m 2 . 

We have mentioned already, that generally speaking one has to consider a total chain of 
compensation equations including connected Green functions with six, eight, etc. legs. Note, 
that corresponding equations will contain inhomogeneous parts, expressed in terms of Green 
functions of lower order, and homogeneous parts, being proportional to the corresponding 
form-factor, e.g. Fq with six legs. Assuming our result the connected four-leg Green function 
be zero, we come to the conclusion, that inhomogeneous part of equation for F& is zero, so 
trivial solution F 6 = inevitably exists. The analogous considerations lead to conclusions 
on possibility of existence of trivial solutions of all higher Green functions. One may, of 
course, study possibilities of existence of nontrivial solutions as well. However, the purpose 
of the present work is to show that even though one nontrivial solution does exist, so we 
rely on following variant: nontrivial solution for four-leg connected Green function and 
trivial solutions for all higher connected Green functions. The consideration of compensation 
equation for Green function with two legs, which defines mass of the scalar field will be 
performed particularly later on. 

The next step of study should include non-linear equation with account of all possible 
diagrams. However this problem evidently do not admit analytic solution. Approximate 
estimate of non-linear corrections to the form-factor's value at zero will be obtained in what 
follows. Maybe future studies will be connected with numerical methods. We are convinced, 
that the experience achieved in finding of the non-trivial solution will help in formulation 
and realization of numerical methods. Presumably result (|37jh which means the existence 
of a solution only for definite relation between dimensional coupling constant and mass of 
scalar field, will be important. 

4 Bethe-Salpeter equation and zero excitation 

It is well-known, that a symmetry breaking is to be accompanied by an appearance of an 
excitation with zero mass PQI21E]- Let us consider this problem in the same approximation. 



12 



While constructing an equation for a bound state one has to keep in mind, that here genuine 
interaction (|13|) acts, that one, which is referred to the interaction Lagrangian and remains, 
of course, not compensated. Bethe-Salpeter equation for a massless bound state of two scalar 
fields in this case has the form 

Gn 3 A' G 2 7r 6 A A' G 2 tt 6 / 1 f x , T/N , 

* {x) = 2(^F " "^F" + r 20^ I v ny)dy + 

3 f x f x 

+ ax~ J y2 ^^ dy + 31ogx J y^(y) d v + 

px px poo 

+ 3x\ogx / ^>{y)dy + A I y^(y)dy + 3 / y log y dy + 



+ x r(4 + 3] ogy )*(y)dy + ^- f°°—dy-^ r^T-dy)- (38) 
J X ^ J X y ^ j x t J 

POO 

A' = / *(y)dy. 
Jo 

Comparing this equation (J3%|) with compensation equation ()19J) . we see the main difference in 
the sign afore the kernel of the integral equation. Remind once more, that the compensation 
equation is the condition of vanishing of the total expansion in G in the modified free 
Lagrangian in expression (JTTj) and therefore terms of the first and of the third orders are 
situated in the same part of equation, e.g. in the left-handed one, whereas in the Bethe- 
Salpeter equation the corresponding terms are situated in different parts of equation. 

The sign afore the kernel is very important. This means, that in a differential equation 
sign afore j3 changes as well 

x 4- + 2 }(x 4- + {) ( x 4-) (x 4-) (x 4- - 1) (x 4- - 1 



d/OC ' \ doc / > d/3y / > doc I \ dioc / > d*oc 

x (4- 2 )(4- 3 )-^ 2 )* w = - (39) 

One easily see, that due to absence of term being proportional to x log x in the inhomoge- 
neous part boundary conditions are the following 

a-2 = o_i = a i = an = . (40) 

The change of sign afore f3 leads to changing of asymptotes at infinity 

V k (x) ~ x~ 3/8 exp^4(/3x 2 ) 1/8 exp(^-^ ; k = 1, 2, 8. (41) 

Now we have three decreasing asymptotes (k = 3, 4, 5), two oscillating ones with power 
decreasing (k = 2, 6), and the remaining three are increasing. Using the first five solutions, 
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which allow a definition of integrals at infinity, we fulfill four boundary conditions at zero (|40|). 
AS a result we obtain the following solution of equation (JHHj) 



where constant A is defined by normalization condition of a Bethe-Salpeter wave function. 
Direct calculation [131 leads to result A' = 0, so the inhomogeneous part of equation ([38)1 
vanishes. Thus we have shown, that the equation for a bound state with zero mass has a 
solution. 

The solution being obtained proves the existence of zero mass excitation P3 |2J El m the 
model. Of course definition of a Bethe-Salpeter equation itself is possible only provided a 
non-trivial solution of a compensation equation to exist and thus interaction (J13|) to act. 
The obligatory correspondence between a non-trivial solution of a compensation equation 
and an existence of a zero excitation thoroughly corresponds to Bogolubov quasi-averages 
approach [2|. 

It is interesting to note, that with taking into account of three-fold interaction g 3 in 
the kernel of equation ()38|) the mass of the bound state becomes non-zero. One easily 
understands this, because interaction (|38|) itself leads to dimensional parameter A 3 being 
present and thus the scale invariance being already broken. 

5 Compensation equation for scalar field mass 

Let us look at interaction Lagrangian (|12|) . The mass term there is quite improper. To 
solve the problem one has to formulate a compensation equation for Green function with 
two scalar legs. Let us consider this equation taking into account solution (}3"2"j) and three- 
fold interaction. The compensation equation means nullification of total contribution of 
interaction (J12j) to the mass. In the first approximation the contribution of the four-fold 
interaction is described by the first order diagram "bubble" and that of the three-fold one 
is represented by simple one-loop diagram (see Fig. 2). 



tt(x) = AGf s {(3x 2 /2 s | 3/2, 1, 1/2, 0, 1/2, 0, -1/2, -1) ; 



(42) 



2 

m 




+ 




Fig. 2. Compensation equation for mass of the scalar field. 
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Putting momenta of the external legs to be zero, we have for "bubble" diagram just 
solution (|32|) in the vertex. As a result we obtain the following compensation equation for 



scalar mass 



G f Ftf)#q g 2 

m = 



\2 



(2tt) 6 J q 2 + m 2 (2tt) 6 J (q 2 + m 2 

C r°° Cm 2 
' ' ydy(F (y) + AF(y)) + ^— , I dy F (y) - 



2(4tt) 3 J » yv uvyy ^" 2(4vr) 

g 2 r y 2 dy 



2(4vr) 3 7 (y + m 



2\2 



(43) 



Here in "bubble" diagram we perform an expansion in m 2 and take into account the zeroth 
and the first orders of the expansion. By direct evaluation with the aid of expressions (|3*Tl 
1341 l33|) we obtain that the zeroth order terms is zero and the first order term is equal to 3 m 2 . 
The loop, which is described by the last term in (J43j) . quadratically diverges. Note, that in 
the initial theory (jSJ) we introduce some cut-off A3, which corresponds to a physical limitation 
of a region of applicability of the theory. As a result we have the following compensation 
equation for the mass provided m -C A 3 



m 2 = 3 m 2 — A3. (44) 



2(4tt) 



Emphasize, that for the trivial solution G = the first term in the right-hand side of 
equation is absent and we have a negative mass squared, i.e. a tachyon solution. For 
the non- trivial solution we have 

m 2 = g2 - A 2 (45) 
4(4vr) 3 3 v ; 

It is well-known, that a scalar tachyon leads to instability for small fields. Therefore the 
restoration of the normal sign of the mass squared, which is achieved provided the non- 
trivial solution is valid, corresponds to a transition to a more stable state. 

So the value of the scalar mass is defined in terms of initial parameters of the theory g 
and A3. The value of parameter Y (p?Tj) gives the relation of the mass and of the coupling 
constant G of the four-fold interaction. Thus all the parameters entering into the non-trivial 
solution are defined in terms of the initial ones. 

Note that the initial cut-off A 3 corresponds to some boundary energy, which provides 
real physical cut-off of the corresponding integrals. In the physical four-dimensional space- 
time it may be for example the Planck energy 1.22 • 10 19 GeV. Note also that in realistic 
models without elementary scalars (see e.g. 0, [Zj) quadratic divergences in mass of the 
elementary fields are absent. One should expect the expressions similar to (|44|) also would 
lead to relations, which connect the theory parameters with a boundary energy (e.g. the 
Planck one), which enters into logarithmically divergent terms. 
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The final result for effective Lagrangian of the theory after the symmetry breaking occurs 
is the following 

— F(xx,X2, x 3 , x 4 ) 0(xi)0(x 2 )0(£ 3 )0(>4) ; (46) 



2 dx^ dx^ 4! 

where form-factor F is the solution of the compensation equation. 

6 Estimate of non-linearity influence 

Till now our results were obtained in the framework of the linear approximation. The 
decrease of the form-factor at infinity indicates an applicability region of the approximation. 
It evidently is incorrect for large momenta variables because the effective coupling constant 
becomes too small in comparison to constant G, which was used to define the kernel of the 
integral equation. We can roughly take into account an influence of a non-linearity, using 
the following procedure. Let equation (|2*TJ) be valid for small x (we put m 2 = 0). 

d 4 f 2 ^ ( 2 rv ^ o F ( x ) o 2G2 

F{x) ))=-^-> P = jAvf- (47) 

We use this equation with the corresponding boundary conditions (|25|) for x < x , whereas 
for x > Xq one has to take into account a non-linearity. Let us draw attention to the fact, 
that (3 is proportional to G 2 i.e. it contains the form-factor squared. Therefore for x > Xq 
instead of (pT7j) we use the following equation 



,r 2 



In this approximation we have correct behaviour of right-hand sides at small (}4Tj) and at very 
large (}4*H|) values of x. In the intermediate region there is a tear in the rhs. at x = xq. This 
means that the eighth derivative tears at this point. As we shall see soon the form-factor 
and its derivatives up to the fifth order have to be continuous. 

Let us introduce variable y = ^/]3 x. One easily sees that for y — ► oo equation (JIB)) 
defines the following decreasing asymptotics 

FM ^7>-ww + rm<? + ' (49) 

where b is a constant. At the same time equation (|47p with account of boundary conditions 
has the following solution in region (0, yo) 

/2 

F(y) = ^G^(y 2 /256|3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) + 
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+ Ci G^(y7256 I 3/2, i, 1/2, 1/2, 0, 0, -1/2, -1) + 

+ C 2 G^(y 2 /256 | 3/2, 1, 0, 1/2, 1/2, 0, -1/2, -1) + (50) 

+ C 3 G^(y 2 /256|3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) + 

+ C 4 G^(y 2 /256 1 1, 3/2, 1/2, 1/2, 0, 0, -1/2, -1) ; 

where Cj are constants. The appearance of the additional terms with these coefficients 
multiplied by Meijer functions increasing at infinity is due to the fact, that now the decrease 
at infinity is provided by asymptotics (jHU) and thus in region (0, yo) we have to use all 
solutions of equation (|47jl . which fulfill the boundary conditions at zero. The first line here 
is solution (|31jh which was obtained earlier. Let us begin a sequential account of the new 
terms starting from the zero approximation, in which in region (0, yo) we have this old 
solution, i.e. all Ci = 0. This solution is matched to solution (}4*9|) in point yo. It will come 
clear, that in expression (f4*9~|) an account of the first term is sufficient. Then from continuity 
of the function and of its first derivative we obtain the following set of equations 

^Oy 7256|3/2, 1, 1/2, 1/2, 0, 0, -1/2, -1) - \ = 0; 

^ GM/256 | 3/2, 1, 1, 1/2, 1/2, 0, -1/2, -1) - \ = ; (51) 
Solution of the set: 

yo = 8.4980; b = 7.5055. (52) 

The second term in asymptotics (}4*9*|) at yo comprises 7.7 ■ 10~ 6 times the first one, that 
justifies the account of the first term only. The value of the form-factor at zero does not 
change F(0) = 1.1107. 

Now let us take into account two additional terms in (|50|) with coefficients C\ and C 2 , 
which for small y give larger contribution than the remaining two terms. In this case we have 
to match values of the function and of its derivatives up to the third order. One obtains the 
set of four equations with aid of rules of differentiation of Meijer functions [12J. Its solution 
reads 

y = 17.635; 6 = 9.410; C x = 0.0166; C 2 = -0.0538 . (53) 
The value of the form-factor at zero becomes the following 

F(0) = — + — = 1.0936. (54) 

4 7T 

Now let us take into account terms with coefficients C3, C4. We consider them and 
deviations from solution (|53|) as well to be small. Then matching the function and its 
derivatives up to the fifth order, we obtain a set of six linear equations leading to the 
following solution 

Ay = 1.457, Ab = 1.032, AC X = -0.0094, 

AC 2 = 0.0223, C 3 = -0.0249, C 4 = 0.0136. (55) 
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Substituting the last result into we have 

F(0) = 1.1007. (56) 

The sequence of numbers 1.1107, 1.0936, 1.1007 for value F(0) demonstrates stability of the 
result in respect to contribution of non-linear corrections 

7 Conclusion 

Grounding on the results being obtained we express a hypothesis, that in the model under 
consideration a nontrivial solution does exist, which breaks the initial scale invariance and 
leads to a spontaneous appearance of effective interaction in Lagrangian ([46]) . acting in 
a restricted region of the momenta space in accordance with the value of parameter G. 
Effective form-factor F(p) decreases exponentially with oscillations for p 2 — > ± oo, i.e. both 
for space-like and time-like momenta. We confirm the existence of a zero mass excitation, 
which has to be present for an occurrence of spontaneous symmetry breaking. 

We start with the asymptotically free theory of a scalar field (in a six-dimensional space), 
and we obtain as a result the definite theory with interaction breaking scale symmetry. New 
dimensional parameters G^ 1 ^ 2 and m are proportional to parameter A3, which defines the 
initial asymptotically free interaction. Let us emphasize once more, that the interaction 
being obtained is an effective one, that first of all is reflected in a presence of form-factor 
F(p), which is just the solution of compensation equation of N.N. Bogolubov method. At 
momentum infinity the theory becomes asymptotically free again. 

It is quite important, that the problem under consideration has a consistent solution only 
provided triple interaction g 3 is acting. Really, albeit compensation equation (fTS]) contains 
no contribution of this interaction, the non-zero scalar field mass appears only for g ^ 0. If 
it is not the case the value of form-factor at zero F(0) is not unity. In general one can not 
exclude a possibility of condition F(0) = 1 being fulfilled for m — 0. However the experience 
obtained in considering the present problem shows that this condition could be fulfilled only 
provided the model has very peculiar properties. As a matter of fact the problem under 
consideration is defined not by compensation equation ([!%]) only, but by set of equations (fTSl 
|4"3*]) . which explicitly contains a contribution of triple interaction g4> 3 . 

It should be noted, that a possibility of a nontrivial solution strongly depends on the 
choice of the theory. This may be demonstrated by comparison of different signatures of the 
six-dimensional space-time. Namely if one instead of signature 1 + 5 will choose signature 
3 + 3, then in definition © of a transition to Euclidean coordinates the sign afore td 6 p 
changes. As a result all signs change for one-loop integrals. For four-fold interaction we 
restore all previous results by simple substitution G — > — G. However the one-loop integral 
with two three-fold vertices inevitably changes sign and relation (|44j) leads to tachyon mass. 
So we come to the conclusion, that for signature 3 + 3 only the trivial solution G = is 
stable. 
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Of course, we base our conclusions only on exact solutions of approximate equations. 
However it is possible, that qualitative properties of solutions, which manifest themselves 
in the model problem, will be quite useful in study of problems of spontaneous symmetry 
breaking in more realistic cases, when there is no hope for analytic solution of corresponding 
equations and what is possible to apply are just numerical methods. Attractive qualitative 
results are the existence of relations between parameters of the problem and the natural 
appearance of small parameter Y (|37|). The essential result is connected also with the 
conclusion on the stability of the non-trivial solution. The estimate of non-linearity contribu- 
tion, which does not lead to decisive change of properties of the solution, provides additional 
argument on behalf of the present approach. 

The resulting theory is non-local and the question might arise, whether the general 
principles of unitarity and causality are here valid. The initial theory (jSJ) quite corresponds 
to these principles. One should expect, that its solutions, non-trivial ones as well, have 
also to fulfill these conditions. Therefore one can consider the present example as a step in 
direction of formulating of a consistent non-local theory. Basing on results of the present 
work we may assume, that such theory can be consistent not for an arbitrary form-factor 
but for the one, which follows from a non-trivial solution of an initially local theory. 

Without any doubt a possibility of spontaneous appearance of an effective interaction, 
containing a dimensional parameter, is of great interest for studies of problems beyond the 
standard theory. In particular, the phenomenon of a spontaneous appearance of an effective 
interaction, provided it to occur in a genuine physical theory, e.g. in the electroweak theory, 
might essentially promote our understanding of bases of the theory. Some considerations 
in connection with this aspects are described in work [HI EI- A subsequent results in this 
direction and their connection with possible deviations from predictions of the standard 
theory will be presented elsewhere. 

The work is supported in part by grant RFBR 01-02-16209 and by grant "Universities of 
Russia" UR.02.03.002. 
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A Appendix 

Here formulas of integration by parts of expressions entering in equation (fT9"|) are presented 

+2 *IM(^>)) - 4(* 2 I^>)) + 

J a (y + mY " dx'\ dx l \ 1 ')) dx 2 \ dx l \ K ')) 
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+ A8x^—(x 2 F(x)) - A8^-(x 2 F(x)) + 48a 
dx 2 V / ax V / 

y 5 F(y) 



d 







o (y + m 2 )' 



dy 



d f n d / > 
-X"-r^\X -r-r\X F{x / 



dx 3 V dx 4 



5x 



d 2 / d 4 / 

4 a_ 2 a_ 2 



20x 3 -^-(x 2 ^—{x 2 F(x" 
dx V ax 4 



d 4 

60x 4 — (x 2 F(x) ) - 120x 3 -f^(x 2 F(x)) + 



dx 4 



dx 2 V dx 4 V 
d 3 



x^f(x; 



dx 3 



+ 360x 2 ^(x 2 F(x)) - 720x^(x 2 F(x)) + 720x 2 F(x) - 720 

m dy „ _! 

(y + m 2 ) 2 dx 3 V dx 4 

yF(y) d 3 , , d 4 , 



a_ 2 



^ (x 2 -^(x 2 F(x) 



x (y + m 2 f 



X dxH X dxA* 



-7^7 TT (^ F(x) 

dx 2 V dx 4 V 



OO 2 



y log !/%) 

(y + m 



dy = x 2 log x -^-g (x 2 -^j-^ ^x 2 F(x)jj — (2xlog x + x) x 



2\2 



^4^^)))+ P log- + 3) 

- 2 4( i!f(i) ) + 4( i!f(s 



(y + m 2 ^ 2 

s(^(^>)); 



2x ^U x2F{x - 



+ 



(y + m 



2\2 



dy = x 3 log x -— (a; 2 ^— j ^x 2 F(x) j j — (3x 2 log x + x 2 ) x 



d 4 / o „, A „ d 3 



a; . 



-(61ogx + ll)x 2 — (x 2 F(x)) + 6x — (x 2 F(; 
For equation ()38|) one has to change sign afore (3 and change F(x) for \l/(x). 



6-^(x 2 F(x) 
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